Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN KTpX style file v2.2) 



Resolving the pulsations of subdwarf B stars: HS 0039+4302, 
HS 0444+0458, and an examination of the group properties of 
resolved pulsators 

M. D. Reed, 1 * D. M. Terndrup, 2 A.-Y. Zhou, 1 
C. T. Unterborn, 2 D. An, 2 and J. R. Eggen 1 

Department of Physics, Astronomy and Materials Science, Missouri State University, 901 S. National, Springfield, MO 65897 USA 
2 Department of Astronomy, The Ohio State University, 140 W. 18th Ave., Columbus, OH 43210 USA 

Accepted Received 

ABSTRACT 

We continue our program of single-site observations of pulsating subdwarf B (sdB) 
stars and present the results of extensive time series photometry of HS 0039+4302 
and HS 0444+0458. Both were observed at MDM Observatory during the fall of 
2005. We extend the number of known frequencies for HS 0039+4302 from 4 to 
14 and discover one additional frequency for HS 0444+0458, bringing the total to 
three. We perform standard tests to search for multiplet structure, measure amplitude 
variations, and examine the frequency density to constrain the mode degree I. 

Including the two stars in this paper, 23 pulsating sdB stars have received follow- 
up observations designed to decipher their pulsation spectra. It is worth an examina- 
tion of what has been detected. We compare and contrast the frequency content in 
terms of richness and range and the amplitudes with regards to variability and diver- 
sity. We use this information to examine observational correlations with the proposed 
k pulsation mechanism as well as alternative theories. 
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1 INTRODUCTION 



Subdwarf B (sdB) stars are horizontal-branch stars with masses ~ 0.5M M , thin (< 10 2 M^) hy 



drogen shells, and temperatures from 22 000 to 40 000 K (|Heber 



1984; 



Saffer et al. 



1994). Pulsat- 



ing sdB stars come in two varieties: short period (90 to 600 seconds), and long period (45 minutes 
to 2 hours). This work concentr ates on the short-pe riod pulsators, which are named EC 14026- 



2647 stars after the prototype ({Kilkenny et al. 



1997b : they are also known as V361 Hya stars or 



sdBV stars. They typically have pulsation amplitudes near 1%, and detailed studies reveal a few 
to dozens of frequencies. The longer period pulsators ar e known as PG 17 16 pulsators after that 
prototype and typically have amplitudes less then 0.1% ( Green etal] 2003 ). They are also cooler 
then the E C 14026-type pulsa tors, though there is some overlap, and they are most likely g— mode 



pulsators ([Fontaine et al. 



2006). 



Asteroseismology of pulsating sdB stars can potentially probe the interior structure and provide 
estimates of total mass, shell mass, luminosity, helium fusion cross sections, and coefficients for 
radiative levitation and gravitational diffusion. To apply the tools of asteroseismology, however, 
it is necessary to resolve the pulsation frequencies. This usually requires extensive photometric 
campaigns, preferably at several sites spaced in longitude to reduce day /night aliasing. Gener- 
ally, discovery surveys have simply identified variables and detected only the highest-amplitude 
pulsations, while multisite campaigns have observed few sdBV stars. 

We have been engaged in a long-term program to resolve poorly-studied sdB pulsators, princi 



pally from single-site data. Th is method has proven useful for several sdBV stars (|Reed et al 



2006 



a, 



2007; 



Zhou et al 



2004, 



2006J). Here, we report the results of our observations of HS 0039+4302 



(hereafter HS 0039) and HS 0444+0458 (hereafter HS 0444). HS 00 39 (B = 15.6) and HS 0444 



0stensen et al 



(2001 



here- 



(B = 16.5) were discovered to be members of the EC 14026 class by 
after 001). Their observations of HS 0039 consisted of three runs of 1, 2, and 3 hours, and they 
obtained three ~ 1.5 hour runs for HS 0444. Even with these short runs, they were able to detect 
four frequencies in HS 0039 and two in HS 0444. 001 also obtained spectra of both stars, from 
which they determined T c g = 33 100 K and 34 500 K, and \ogg = 6.0 and 6.1 (with g in cgs units 
of cm ■ s~ 2 ) for HS 0039 and HS 0444, respectively, and that neither star is a spectroscopic binary. 
In §2 we describe our new observations of these stars, in §3 we analyze the pulsation frequencies, 
and in §4 we discuss our findings and apply asteroseismic tests. 

With the addition of these two stars, 23 sdBV stars have received follow-up observations, 
including 18 stars for which our program has directly contributed data. In §5, we will discuss 
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the observational properties of all 23 stars, concentrating on pulsation stability, amplitudes, and a 
comparison with known driving mechanisms. 



2 OBSERVATIONS 

Data were obtained at MDM Observatory's 1.3 m telescope using an Apogee Alta U47+ CCD 
camera. MDM Observatory is located on the southwest ridge of Kitt Peak, Arizona and is operated 
by a consortium of five universities, including the Ohio State University. Images were transferred 
via USB2.0 for high-speed readout; our binned (2 x 2) images had an average dead-time of one 
second. The observations used a red cut-off filter (BG38), so the effective bandpass covers the B 
and V filters and is essentially that of a blue-sensitive photomultiplier tube. Such a setup allows 
us to maximize light throughput while maintaining compatibility with observations obtained with 
photomultipliers. Tables Q] and [2] provide the details of our observations including date, start time, 
run length, and integration time. The observations total nearly 100 hours of data for HS 0039 and 
more than 60 hours for HS 0444. 

Standard image reduction procedures, including bias subtraction, dark current and flat field 
correction, were followed using IRAFLj packages. Intensities were extracted using aperture pho- 
tometry. Extinction and cloud corrections were obtained from the normalized intensities of sev- 
eral field stars. Because sdB stars are substantially hotter than typical field stars, differential light 
curves are not flat due to atmospheric reddening. A low-order polynomial was fit to remove nightly 
trends from the data. Finally, the lightcurves were normalized by their average flux and centered 
around zero so the reported differential intensities are AI = (//(/)) — 1. Amplitudes are given as 
milli-modulation amplitudes (mma), with 10 mma corresponding to 1.0% or 9.2 millimagnitudes. 
Sample lightcurves are shown in Fig. [TJ 



3 PULSATION ANALYSIS 

HS 0039: A quick analysis during observations alerted us that the amplitudes of HS 0039 were 
not stable. This is easily seen in the final nightly reductions, twelve of which are shown in Fig. [2] 
In this figure, we show the pulsation spectra (Fourier transforms; FTs) from adjacent nights, ex- 
cept for the short night on 27 November. Only the frequency near 4270 fiRz appears stable in 
amplitude, while those near 5175, 5482, and 7348 fiRz show substantial variation. We therefore 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Figure 1. Representative lightcurves of HS 0039 (top pair) and HS 0444 (bottom pair). Each panel is one hour in length. The line shows our fit to 
the data. Note that the vertical scale differs for the two plots. 



examined the amplitudes and phases for indications of closely spaced multiplets, which would 
prod uce roughly sinusoidal amplitude v ariations and phase changes near the median amplitude 
(see LDaszyhska-Daszkiewicz et al.l 120051) . Figure [3] shows our non-linear least-squares fits for the 
four dominant frequencies. 

The frequency near 4270 fiRz is the most stable, both in amplitude and phase. The remaining 
three frequencies show significant amplitude variation, but only a little variation in phase. The 
phase for the peak near 5175 fiHz shows an unusual bimodality at the beginning of the campaign, 
with a steady, intermediate value at the end. We therefore separated the data into two subsets 
composed of data from the first seven runs (15 - 21 Nov.) and the last six runs (25 - 30 Nov.). 
Figure H] shows the region near 5175 fiRz for all the November data and the two subsets. The 
FT of the first seven runs allowed us to interpret the phase information: the frequency at 5175 
/iHz is composed of a close doublet separated by ~ 17/iHz. The separation between frequencies 
is just shorter than 1 day (11.6 /iHz) and so nightly runs will not resolve these into two separate 
frequencies. Evidently the timing was just right near the beginning of the campaign that the phase 
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Table 1. Observations of HS 0039 



Run UT Start UT Date Length Integration 

(h:m:s) 2005 (hr) (s) 



mdml 11505 


01:46:00 


15 


Nov. 


7.9 


15 


mdml 11605 


01:17:00 


16 


Nov. 


8.4 


12 


mdml 11705 


01:18:00 


17 


Nov. 


8.4 


12 


mdml 11 805 


01:16:00 


18 


Nov. 


8.5 


12 


mdml 11905 


01:15:00 


19 


Nov. 


8.5 


12 


mdml 12005 


01:23:00 


20 


Nov. 


8.3 


12 


mdml 12105 


01:23:00 


21 


Nov. 


8.2 


10 


mdml 12205 


01:11:00 


22 


Nov. 


8.6 


10 


mdml 12505 


02:35:00 


25 


Nov. 


6.5 


12 


mdml 12605 


01:06:00 


26 


Nov. 


4.0 


12 


mdml 12705 


01:11:00 


27 


Nov. 


0.6 


12 


mdml 12805 


01:11:00 


28 


Nov. 


6.2 


10 


mdml 12905 


01:06:00 


29 


Nov. 


6.4 


10 


mdml 13005 


01:90:00 


30 


Nov. 


3.4 


10 


mdml20905 


01:24:00 


09 


Dec. 


4.7 


8 


mdml21405 


01:14:00 


14 


Dec. 


1.1 


8 



Table 2. Observations of HS 0444 



Run UT Start UT Date Length Integration 

(h:m:s) 2005 (hr) (s) 



hs04mdml 11505 09:54:00 

hs04mdml 11605 09:52:00 

hs04mdml 11705 09:53:00 

hs04mdml 11805 09:52:00 

hs04mdml 11905 09:53:00 

hs04mdml 12005 09:57:00 

hs04mdml 12105 09:48:00 

hs04mdml 12205 10:52:00 

hs04mdml 12605 05:26:00 

hs04mdml 12705 08:59:00 

hs04mdmll2805 07:23:10 

hs04mdml 12905 07:32:00 

hs04mdml 13005 04:26:00 

hs04mdml20905 06:26:40 

hs04mdml21005 01:30:30 

hs04mdml21405 05:45:10 



15 Nov. 


3.0 


15 


16 Nov. 


3.1 


12 


17 Nov. 


3.0 


12 


18 Nov. 


2.9 


12 


19 Nov. 


2.8 


12 


20 Nov. 


2.7 


12 


21 Nov. 


2.8 


12 


22 Nov. 


1.6 


12 


26 Nov. 


6.8 


12 


27 Nov. 


3.1 


12 


28 Nov. 


4.8 


10 


29 Nov. 


4.5 


10 


30 Nov. 


4.7 


10 


09 Dec. 


2.6 


15 


10 Dec. 


9.5 


10 


14 Dec. 


5.3 


12 



switched between the two frequencies of the doublet on alternate nights. This was not the case 
later in the run, although the very last phase might show that the pattern was reestablishing itself. 

The frequency near 5482 /xHz shows a similar variation in amplitude to that near 5 175 /xHz, but 
not the bimodal phase. We grouped the data into the same subsets as in Fig. |U but did not see any 
clear sign of a close doublet. If the amplitude variations were intrinsic to that peak, the discovery 
of frequencie s by prewhitening would be more complicated. We produced an accurate window 



function as in 



Reed et al. 



(|200 6a); the window function is the FT of a noise-free sinusoidal single 
frequency sampled at the same times as the data. The window function matched the multi-peaked 
structure of the temporal spectra around 5482 /jHz, leading us to conclude that the amplitude and 
phase variations are intrinsic to that frequency and it is not a closely-spaced multiplet. We also 
examined the region near 7348 /iRz in the complete data and in the subsets. There is no indication 
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Figure 2. Pulsation spectra of temporally adjacent data runs for HS 0039 plotted at the same scale. 



of closely spaced multiplets, nor in the the phases (Fig. [3]). We conclude that this frequency also 
has intrinsic amplitude variations, but is a temporally resolved frequency. 

We do not see any such indications in the remaining frequencies, though most have amplitudes 
that are too small to be detected in individual runs. Armed with five resolved frequencies, we 
continued simultaneously least-squares fitting and removing peaks (prewhitening) in the combined 
data set. This process is shown in Fig. [51 where the top panel is the original FT and the next three 
panels show residuals after prewhitening by 5, 10, and 14 frequencies (from top to bottom). The 
solid (blue in the on-line version) line in the figure indicates the 4cr noise limit, below which we 
do not fit any peaks. After fitting 14 frequencies, we concluded that any remaining power in the 
residuals was caused by amplitude and/or phase variations and so does not represent remaining 
unfit frequencies. The solution to our fit, listing frequencies, periods, and average amplitudes is 
provided in Tabled 

HS 0444: Early in our campaign, we concentrated on HS 0039, leaving HS 0444 as a secondary 
target. Even from the relatively short runs (Table [2]), we could tell that HS 0444 was a simple 
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Figure 3. Amplitudes and phases of the four largest amplitude frequencies of HS 0039. Frequencies are provided (in /jHz) in each panel. 



Table 3. Frequencies, periods, and average amplitudes for HS 0039. Formal least-squares errors are in parentheses. 



ID 



fl 

n 

/3 
/4 
/5 
/6 
/7 
/8 
/9 
/10 

/H 
/12 
/13 
/14 



Frequency 
(/.Hz) 



Period 

(s) 



4135.559 
4271.481 
4952.858 
5130.474 
5175.466 
5192.660 
5482.319 
5550.170 
5705.677 
5751.316 
5756.651 
7348.444 
7449.256 
7459.765 



(0.042) 
(0.008) 
(0.078) 
(0.036) 
(0.009) 
(0.010) 
(0.016) 
(0.047) 
(0.044) 
(0.081) 
(0.068) 
(0.031) 
(0.064) 
(0.063) 



241.8052 
234.1108 
201.9035 
194.9137 
193.2193 
192.5794 
182.4045 
180.1746 
175.2640 
173.8732 
173.7120 
136.0832 
134.2415 
134.0524 



(0.0024) 
(0.0004) 
(0.0031) 
(0.0013) 
(0.0003) 
(0.0003) 
(0.0005) 
(0.0015) 
(0.0013) 
(0.0024) 
(0.0020) 
(0.0006) 
(0.0011) 
(0.0011) 



Amplitude 
(mma) 



0.85 (0.09) 
4.35 (0.09) 
0.46 (0.09) 
1.00 (0.09) 
3.99 (0.09) 
3.73 (0.09) 
2.16(0.09) 
0.77 (0.09) 
0.82 (0.09) 
0.45 (0.09) 
0.53 (0.09) 
1.15 (0.09) 
0.65 (0.09) 
0.67 (0.09) 



pulsator with just a few frequencies. In December, when HS 0444 was better placed in the sky, 
we obtained longer individual runs to ensure that we did not miss any low amplitude peaks. The 
temporal spectrum of HS 0444 is shown in Fig.[6l The top panel shows the original FT, the bottom 
panel the residuals after prewhitening by three frequencies, and the inset is the window function. 
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Figure 4. Close-up of the region near 5175 /iHz showing prewhitening of the frequency doublet in HS 0039. From top to bottom the panels are the 
window functions, the original FTs, and the FTs prewhitened by one, and two frequencies. The left panels are the complete data, the middle panels 
are the first seven data runs combined, and the right panels are the last six data runs combined. 



The frequencies, periods and amplitudes from our least-squares solution are provided in Table 
HI In Fig. |71 we show the amplitudes and phases of the two frequencies which were detectable 
each night. The amplitudes and phases are relatively stable and do not show any indications of 
additional unresolved frequencies. 
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Figure 5. Temporal spectrum for the complete HS 0039 data set showing the prewhitening sequence. The top panel shows the original FT while 
the remaining panels show the data prewhitened by 5, 10, and 14 frequencies, respectively. The prewhitened frequencies are indicated by arrows 
and the solid (blue in the on-line version) lines indicate the 4<r detection level while the inset is the window function. 



Table 4. Periods, frequencies, and amplitudes for HS 0444. Formal least-squares errors are in parentheses. 



ID 


Frequency 


Period 


Amplitude 




(fiHz) 


(s) 


(mma) 


/I 


5902.511 (0.012) 


169.41940 (0.00035) 


2.5 (0.2) 


/2 


6553.484 (0.037) 


152.59056 (0.00086) 


0.8 (0.2) 


/3 


7311.728 (0.003) 


136.76657 (0.00005) 


11.1 (0.2) 



4 DISCUSSION 

4.1 Comparison with the discovery data 

The goal of our observational program is to resolve the pulsation frequencies for asteroseismic 
analysis. For the sake of comparison with the discovery data, we ca lculate the temporal r esolution 
as 1/At with At being the extent of the observations in time (see iKilkenny et al.ll 1999T) . For the 
discovery data, we can determine the temporal resolution from information provided in 001, and 
estimate the detection limit as twice the top of the noise level in their FTs outside of the pulsations 
and window functions. 
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Figure 6. Temporal spectra for HS 0444. The top panel shows the original FT, the bottom the residuals after prewhitening by three frequencies, 
and the inset is the window function. The solid (blue) line is the 4<r detection limit. 

For HS 0039, our observations have a temporal resolution of 0.8 /iHz (excluding the De- 
cember runs), which is « 6.5 x better than the discovery data. The detection limit is 0.4 mma, 
which is about 5x better than in 001. It is difficult to determine whether the frequencies in 
HS 0039 have changed since the discovery observations. 001 's reported frequencies (5.14, 5.48, 
and 4.27 mHz) roughly correspond to our /4 or /5, /7, and f2, respectively, but 001 's last fre- 
quency of 5.21 mHz, detected in only a single observing run, is not found in our data. As this 
frequency is wedged between higher-amplitude frequencies, and 001 's window function is com- 
plex, it is difficult to judge the significance of its detection, though it does look reasonable in their 
figures. However, the main difficulty is the abundance of frequencies, many of which have low 
amplitudes. The discovery data only detected four frequencies, whereas we detect 14. We estimate 
001 's noise limit as ps 2 mma; if our limit were as high, we also would have detected four fre- 
quencies. It is therefore not possible to ascertain the long-term stability of the pulsations from just 
these two sets of data. 

Additional data on HS 0039 were obtained using the ULTRACAM multicolor instrument on 
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Figure 7. Amplitudes and phases for frequencies / 1 and /3 of HS 0444. 



the 4.2 m William Herschel Telescope in 2002 by Jeffery et al. (2004; hereafter J04). Two long 
(w 8 hr) runs with very high signal-to-noise were obtained on consecutive nights. We recover 
J04's eight frequency detections, six of which are independent. Two are aliases of f7 and one is 
an alias of /13. J04 used multicolour photometry to estimate that /8 is an £ = 4 mode. Though 
geometric cancellation should reduce amplitudes for higher degree modes, the (1, 771 ) = (4, 1) 
mode can be reduced by as little as 30%, depending on orientation (IReed et alj 120051) . Since /8 
is 18% the amplitude of f2 (the highest amplitude frequency), there is no problem in ascribing 
it as an i = 4 mode. Additionally, as this region of the FT is relatively uncrowded, the data 
obtained by J04 should have been sufficient to make this determination. Although they do not 
claim this identification with certainty, it seems reasonable and would certainly be worth additional 
multicolour photometry or time-series spectroscopy for confirmation. 

For HS 0444, our observations have a temporal resolution of 0.4 /iHz, which is w 14 x better 
than the discovery data. The detection limit is 0.6 mma, which is about 3x better than in 001. 
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Additionally, we have recovered the two frequencies detected in the discovery data, to within the 
errors, and uncover a single new frequency, at an amplitude below their detection limit. 



4.2 Constraints on the pulsation modes 



In addition to improving the known pulsation spectra of these stars, we wish to place observa- 
tional constraints on the pulsation modes. The modes are mathematically described by spherical 
harmonics with three quantum numbers, n (or k), £, and m. Rotation can break the m degeneracy 



usefu l tool for observatio nally constraining pulsation degree (see 



2004 



so multiplet structure is a very 


Winget et al. 


1991 


•> 


Reed et al. 



OToole et al.l l2004). 



For slow rotators, like most sdB stars are thought to be (|Heber. Reid. Werner 



1999 



2000), 



rotationally-split multiplets should be nearly equally spaced in frequency. Such structure is, how- 
ever, seldom observed in sdBV stars and in neither case are multiplets detected in our observations. 
Even though HS 0039 has 14 frequencies, the frequency spacings are not regular. Instead, the spac- 
ings are distributed from 5 to 3325 /iHz, with no obvious groupings. For HS 0444, there are only 
three detected frequencies, but the spacings are not similar, so there is no multiplet structure in this 
star. 

Another tool that can be used is the density of frequencies within a given range. In resolved 
sdBV stars, we sometimes observe many more pulsation modes than £ = 0, 1, and 2 can provide, 
independent of the number of inferred m^0 frequencies. Higher £ modes may be needed, but if 
so they must have a larger amplitude than is measur ed because of the large degree of geometric 



cancellation (ICharpinet et al. 



per i degree per 1000 /xHz (ICharpinet et al. 



2005 



Reed et al.||2 005). A general guideline would be one n order 
( 2002() find an average spacing near 1440 /iHz), so 



the temporal spectrum can accommodate three frequencies per 1000 /uHz without the necessity of 
invoking high-£ values if no multiplet structure is observed. Filling all possible m values, the limit 
becomes nine frequencies per 1000 /^Hz.. 

Obviously there is no need to invoke high-£ modes for HS 0044. Between 4900 and 5800 /iHz, 
HS 0039 has 9 of its 14 frequencies. Therefore, HS 0039 has too large a frequency density to 
exclude £ ^ 3 modes, particularly with the absence of any obvious multiplets. This supports the 
identification of an £ = 4 mode by J04. 
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5 GROUP PROPERTIES 
5.1 Data sources 



1997), 



Since the discovery of the EC 14026 class of pulsating sdB stars in 1997 (jKilkenny et al 
there have been three areas of emphasis for observations: 1) to discover more pulsators; 2) to re- 
solve the pulsations using long time-base campaigns, sometimes at multiple sites; and 3) to obtain 
high signal-to-noise observations over short time intervals. In recent years, multicolour photome- 
try and time-series spectroscopy have been obtained as additional tools for mode identification. 

For this paper, we will concentrate on the second point above, and examine pulsators for which 
a considerable effort has been expended to resolve the pulsation frequencies. We do this because it 
has been our area of emphasis, we feel that it is an important component in applying asteroseismol- 
ogy to sdB stars, and most importantly, we have data for all the stars except those from Kilkenny 
et al. (2002, 2006a, 2006b). Though this will not be a complete sample of sdBV follow-up obser- 
vations, we can perform uniform tests upon them (except as noted above) for intercomparison. 

Table [5] provides a list of studies that have thoroughly investigated the pulsations of sdBV 
stars. Some stars, such as Feige 48, PG 0014, PG 1219, and PG 1605, have received extensive 
observations over the course of many years, while most have only been observed during a single 
campaign. Column 1 of that table lists the full name of each pulsator. Columns 2 and 3 display 
the dates of observations and the number of hours observed. The final two columns display the 
observing sites and references for each star. The references in Table|5]are the basis for our analysis 
below, but are not necessarily a complete record of the detailed observations of each star. 



5.2 The pulsation content 

In Table [6j we assemble data on the various sdBV stars. Column 1 lists an abbreviated name, 
which we will use hereafter. Columns 2-4 give the total number of detected frequencies, and 
the numbers with high and low amplitudes (these and other quantities in the table are discussed 
more thoroughly below). In columns 5 and 6, we show the temporal resolution 1/At and the 
noise limit in mma as inferred from the analysis in each paper (§4.1). Column 7 displays our 
judgment about whether the frequencies were completely resolved. Of the 23 stars in Table [71 18 
are likely resolved. The other five either have too little data (EC 20338 and PB 8783), or have 
(some) pulsation amplitudes that are variable on time-scales too short for frequency resolution 
(BA09, KPD 1930, and PG 0048). The effective temperatures and gravities are listed in Columns 
8 and 9; quantities in parentheses are estimates (described in §5.3). The next two columns display 
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Table 5. List of follow-up observations of pulsating sdB stars. A dagger+ indicates observations that we were involved in and NA indicates 
information that was not available. Observing sites: 1) Suhora 0.6 m; 2) Baker 0.4 m; 3) CTIO 1.5 m; 4) SAAO 1.9 m; 5) Fick 0.6 m; 6) MDM 
1.3 m; 7) McDonald 2.1 m; 8) McDonald 0.9 m; 9) MDM 2.4 m; 10) Whole Earth Telescope Campaign, 11) Multisite campaign by Reed et al., 
12) Other campaign. Column 5 references are: a) Baran et al. 2005; b) Reed et al. 2006b; c) Kilkenny et al. 2006a; d) Kilkenny et al. 2006b; e) 
Reed et al. 2004; f) this work; g) Reed et al. 2006a; h) Silvotti et al. 2002a; i) Reed et al. 2006c; j) Zhou et al. 2006; k) O'Donoghue et al. 1998a; 1) 
Vuckovic et al. 2006; m) Reed et al. 2007; n) Kilkenny et al. 2002; o) Harms, Reed, O'Toole 2006; p) Silvotti et al. 2006; q) Kilkenny et al. 2003; 
r) Kilkenny et al. 1999; U) unpublished. 
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the total power in the resolved frequencies and the largest detected amplitude. The last column 
provides references for the spectroscopically determined values of log g and T eff . 

Figure [8] displays log g against T eS for the 20 stars in Table [6] with previously determined 
values. These are shown as filled circles. In addition, PG 1716-type pulsators are shown as filled 
(blue) triangles, and non-pulsators are shown as open circles. The solid (black) line is the zero-age 
helium main sequence with masses marked, the dash ed line is the zero-age extended horizontal 



branch, and evolutionary tracks (from 



Reed et al. 



20041) are shown as solid (blue) lines. The coolest 



track has a hydrogen envelope thick enough for shell fusion, while the hotter two do not. 

Figures|9]and[10]show schematic representations of the temporal spectra of all 23 stars, ordered 
by log g, with the three stars without spectroscopically constrained values at the end. The dotted 
arrows (blue in the electronic version) indicate frequencies which are only observed occasionally. 
PG 0014 has dashed arrows (green in the electronic version), indicating frequencies that were 
only observed using ULTRACAM. To make low-amplitude frequencies visible in the plots, the 
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Figure 8. A log g — T c ff diagram showing the location of sdB stars. Filled circles are EC 14026-type pulsators, filled (blue) triangles are PG 1716- 
type pulsators, and open circles are non-pulsating sdB stars. The dashed line is the zero-age extended horizontal branch, the black line is the zero- 
age helium main sequence with open circles indicating the total mass. The evolutionary tracks (solid blue lines) were produced using ISUEVOS 
jReed et alj2004) . 

vertical axes may begin below zero and are scaled so that the highest amplitude peak touches the 
top line, except for PG 1605, BA09, PG 1325, and EC 20338. All of the latter have one high- 
amplitude peak that would make the others too small on the plot. Those frequencies are indicated 
with arrows which pass beyond the top of the plot, though still not to scale. Note that all panels 
are plotted at different scales. BA09 shows a mixture of short- (EC 14026-type), and long-period 
(PG 1716-type) oscillations; only the former are shown. 

These two figures show the enormous variety of amplitudes and frequencies detected in sdBV 
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Table 6. Pulsation properties of EC 14026-type pulsators for which follow-up data has been obtained. References for the spectroscopic measure- 
ments are as follows: 1) Oreiro et al. 2004; 2) Koen et al. 1999a; 3) O'Donoghue et al. 1997; 4) Koen et al. 1998b; 5) 0stensen et al. 2001a; 6) 
0stensen et al. 2000b; 7) Billeres et al. 2000; 8) Koen 1998; 9) Brassard et al. 2001; 10) O'Donoghue et al. 1998b; 11) Koen et al. 1999b; 12) Sil- 
votti et al. 2002b; 13) Kilkenney et al. 1998; 14) Koen et al. 1998a; 15) Silvotti et al. 2000; U) Unpublished. Parameters in parentheses are inferred 



using color and pulsation frequencies. NA indicates information that is not available. BA09 also has longer period (PG 1716-type) pi 
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stars. There are only four high- amplitude (here A > 20 mma) pulsators known and they have a 
great range (for sdB stars) of gravities and temperatures. There are pulsators with 20+ frequencies 
that have similar temperatures and gravities to stars with 5 frequencies (e.g., BA09 & HS 2201). 
If one looks at only the range of gravities from log g = 5.6 to 5.7 (about la in error), there 
are two stars with more than 20 frequencies yet one star with only three frequencies (but higher 
amplitudes!). 

5.3 Observational tests and trends 

In an attempt to bring order to the class as a whole and to find trends in the observational prop- 
erties, we have organized the data in several ways which have benefited studies of other variable 
stars. In this subsection we will show the results along with some motivation, but leave in-depth 
interpretations to the next subsection. 

Frequency groupings: Our first arrangement was to put the frequencies shown in Figs. |9] and 
[TOl onto a common frequency scale and make a correction for gravity. Because p— mode periods 
are inversely proportional to the square root of the density, we expect that stars with lower log g 
have longer periods (shorter frequencies). This is largely observed in the left panel of Fig. [TT] 
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Figure 9. Schematic representation of the temporal spectra of sdBV stars with follow-up data ordered by log g (from lowest to highest). Blue 
arrows indicate frequencies that are not regularly detected and green arrows for PG 0014 indicate frequencies only detecting using ULTRACAM. 
Shown here are the stars with lower gravities. 



In the right panel, we have rescaled pulsation frequencies by l/g 75 (assuming constant mass) 
which is an adjustment for both size and density using just log g. For the three sdBV stars without 
measured gravities, we estimated the gravity from the position of the shortest frequency compared 
to pulsators with known gravities. These assumed gravities are provided in parentheses in Fig. ITTI 
and Table [6] Such a study of pulsati ng white dwarf s revealed groups of frequencies which could 



then be related to individual modes (IClemens 



1994). However, as evidenced by the summation of 



the right panel, no such groupings occur. The horizontal line just above the summation frequencies 
(right panel of Fig. [H]) shows the effect of an error of log g = 0.05 in the rescaling. It is therefore 
possible that any groups are being smeared out by measurement errors in logg. We attempted to 
correct for this by fixing the lowest modified frequency to a given value, but no corrections or fixed 
reference values of log g show reasonably- separated grouping that could be of use. 

Relative pulsation amplitudes: The line lengths in Figs.l9lthrough[TT1indicate another feature 
that is observed about half the time; one or two amplitudes are significantly higher than the rest. 
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Figure 10. Same as Fig. [9] but for stars with higher gravities. The last three panels are for the stars without a previously determined value of log g. 
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To parameterize this phenomenon, we have denoted frequencies whose amplitudes are within a 
factor of five of the highest amplitude as "high" amplitudes; the remainder are "low" amplitudes. 
The numbers that fall into each category, along with the highest amplitude observed (A max ) for 
each star, are provided in Table [6J The factor of five was chosen so that all amplitudes greater 
than 10 mma for PG 1605 and BA09 would fall into the high category. Additionally, we use the 
published (typically average) amplitudes, even though some vary by large amounts (discussed 
below). 

The left panel of Fig. [12] shows the ratio of high (H) to total (T) number of frequencies against 
the total number of frequencies; the values of which are also provided in Table [6l One might 
suspect that if a star has one high amplitude, all amplitudes are relatively higher and thus easier to 
detect. This appears not to be the case, since the points make a scatter diagram. In the right panel 
of Fig. d2l we show the H/T ratio against log p. It might be expected that gravity plays a factor, 
in that it is easier to pulsate radially at lower gravities. The values of H/T are indeed positively 
correlated with gravity, but the correlation coefficient is only 0.59. 
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Figure 11. Schematic representations of pulsation frequencies organized by log g. The lines used to indicate frequencies have two lengths, with the 
shorter lengths indicating frequencies with amplitudes a factor of 5 smaller than the highest amplitude. The right panel has a frequency correction 
for log g and both panels include a summation of all frequencies at the bottom. 
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Figure 12. Comparison of number of "high" amplitude frequencies (H) to the total number of frequencies (T). Double sized points indicate two 
stars with the same quantities. The right panel compares this ratio with the surface gravity. 
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A more quantitative approach is shown in Fig. [131 where we analyse the distribution of ampli- 
tudes and power in each star. For each star, the frequencies were sorted by increasing amplitude 
and the cumulative distribution of amplitudes (or power) was computed at intervals of 10% of the 
total. This process is shown in Fig.[l4]for the star HS 0039. The 14 frequencies from Table[3]are 
shown as squares with their amplitudes given by the right-hand Y-axis. The circles indicate the 
cumulative fractional amplitude for each frequency, connected by the solid line, with their scale 
given by the left-hand Y-axis. The points for HS 0039 in the left panel of Fig. [T3l correspond to 
the locations in Fig. [141 where the dotted lines intersect the solid line. The top solid line in Fig. IT3l 
connects the points for the 90 th percentile, the next line for the 80 th percentile, and so on. The left 
panel of Fig. [13] displays these points, or amplitude distribution, for the 20 stars with more than 
five frequencies. Their designations are provided along the bottom in order of increasing gravity. 
Low values indicate that the top one or two amplitudes contain a large fraction of the total am- 
plitude (Fig. [14] would show a sharply peaked line). The right panel displays the distribution of 
power which we define as the sum of the squares of the amplitudes; the power and the maximum 
amplitude A max are provided in Tabled On average, the lowest-amplitude half of the frequencies 
contribute about 25% of the total amplitude and 9.6% of the total power. 

If the amplitudes or power were evenly distributed (equal-height amplitudes) then / fr would 



equal Aij 'At (or yPi/Pr)- The corresponding distribution in Fig. [14] would be a straight diagonal 
line. There is a trend in Fig. [13] in that the contours are closest at low gravities and most diffuse 
at high gravities. However, these trends are dominated by four stars (two at each end) and are not 
representative of the majority of the class. In the middle of each panel, there is a large dispersion 
in values, with neighbouring stars transmitting most of their power through one frequency, or 
distributing it more equally. 

In each panel, the dashed line indicates the fractional amplitude or power emitted by the fre- 
quency with the highest amplitude. For PG 1325, the one highest amplitude has 99% of the total, 
while for PG 0048 it is only ~ 12%. The dashed line in Fig. d indicates that for HS 0039, the 
highest amplitude frequency has 20% of the combined amplitudes. While it is difficult to assign 
significance to this plot as the number of pulsators is still relatively low, it is suggestive that in 
lower gravity stars the pulsation energy is channeled into relatively few (or one) frequencies even 
though many frequencies are available. 

Multiplet structure: As the prototype for multiplet structure in pulsating white dwarfs, PG 1 159 
showed that observa tional determination of the modes can lead to tight constraints on the models 



(|Winget et al 



19911) . With detailed studies of 23 pulsating sdBV stars, it would be hoped that a 
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Figure 13. Distribution of amplitudes (left panel) and power (right panel) for pulsators with at least five frequencies. The stellar designations are 
provided at the bottom and are organized by increasing gravity similar. Each colored line indicates a decreasing fraction of pulsation frequencies 
(from top to bottom), with the first four lines labeled. The dashed line indicates the fractional amplitude or power of the highest amplitude frequency 
in each star. 



similar star of this class would have been detected. However, this has not been the case: multiplet 
structure has been conspicuously absent, even from rich pulsators. The only confirmed case of 



multiplets caused by rotational sp litting is Feige 48 (IReed et al. 



2004; 



O' Toole et al 



20041) . while 



another likely candidate is BA09 (IBaran et al.ll2005|) . There are also marginal cases for rotationally 
induced multiplet structure in HS 2201, PB 8783, PG 0014, PG 1047, and PG 1605 none of which 
have been confirmed by additional observations. KPD 1930 and PG 1336 are both known to be in 
short period binaries and frequency splittings are commensurate with the binary period. An initial 
interpretatio n did not attribute these to rotational splitting, but rather to tidal effects induced by the 



companion (IReed et al. 



2006c 



d) 



A possible reason for the lack of observed multiplet structure was proposed by 



Kawaler & Hostler 



(|2005|) . Their picture invokes sharp differential rotation in rapidly spinning cores to modify the fre- 
quency spacings. Exceptions would be for those stars in close binaries where rotations are tidally 
locked. 

Another multiplet pattern that has emerged lately is the "Kawaler-scheme" (Kawaler et al. 
2006) a purely mathematical formalism based on an asymptotic-like relationship for the frequen- 
cies: f(i,j) = f + ix 5+j x A where i has integer values, j is limited to values of — 1 , 0, and 1, 5 
is usually a small spacing and A is usually a large spacing. Their Table 3 indicates the significance 
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Figure 14. Distribution of amplitudes for HS 0039. The frequencies are organized from lowest to highest amplitude with the amplitudes indicated 
by filled squares (on the scale of the right-hand axis). The cumulative fractional amplitudes are indicated by filled circles and connected by a solid 
line. The vertical dotted lines are locate at 10% intervals to indicate the points in Fig. ll3l and the horizontal dashed line (on the scale of the left-hand 
axis) indicates the fractional amplitude of the highest amplitude frequency. 

of their predicted frequencies to those observed in several stars. However, it is also known that the 
Kawaler-scheme does not fit several pulsators (including HS 0039) and so we merely make note 
of the scheme but await a full report in a forthcoming paper. 

Frequency density: As discussed in §4.2, another tool at our disposal is the mode density. 
Although the mode density does not help to assign modes to individual frequencies, we can set 
limits on the number of degrees (£) per order (n) required to create the observed frequency density 



from currently available models of iCharpinet et al 



2004) 
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J2001L hereafter CFB01) and of Reed et al. 



Figure \\5\ shows the mode density, plotted against log g, with the dotted line indicating 3 
frequencies per 1000 fiRz (all i ^ 2, m = modes) and the dashed line indicating 9 frequencies 
per 1000 /iHz (all possible 2£ + 1 m modes). Starred points indicate those pulsators for which 
we have inferred the log g values. Less than 20% of sdBV stars have mode densities below 3 per 
1000 /iHz, while more than 25% have mode densities too high to be reconciled with i ^ 2 modes 
even if all possible m values are used. It therefore seems reasonable to conclude that £ ^ 3 modes 
must be excited. 

Amplitude variability: A criterion outlined in Christensen-Dalsgaard et al (2001; hereafter 
JCD01) and previously applied to several sdBV stars (Pereira & Lopes 2005; Zhou et al. 2006; 
Reed et al. 2006a, 2007) is to compare the average amplitude (A) to the standard deviation of the 
amplitudes, a (A). For stochastically excited pulsations, this ratio should have a value near 0.5. For 
all resolvable frequencies in our target stars, we have calculated both parameters and their ratios 
which are given in columns 3 to 5 of Table [7] and plotted in Fig. [16] Columns 6 and 7 of Table [7] 
list the maximum and minimum observed amplitudes and column 8 gives the maximum time-scale 
over which the observations are sensitive to amplitude variations. In the Figure, triangles (squares) 
indicate frequencies known to have stable (not-stable) phases, circles indicate frequencies with 
ambiguous or no phase information and stars indicate the PG 0048 frequencies, which are known 
to have stochastic-like properties (Reed et al. 2007). Stochastic oscillations do not have stable 
pulsation phases and so frequencies with stable phases should be driven rather than stochastically 
excited. Amplitudes for HS 2201 came from Silvotti et al. (2002a) and no errors were published, 
so no errorbars appear in the Figure. 

Just like their average amplitudes and frequency density, the only conclusion we can draw 
from Fig. [16] and Table [7] is sdBV stars show a large variety of amplitude variations. Of course 
it is known that many classes of pulsators show large amplitude differences and that correlations 
between excitation rates and amplitudes are weak at best, but these results indicate that there is 
no clear separation in amplitude variability for phase-stable and unstable frequencies. As such, it 
follows that the JCD01 criterion is likely not applicable for these stars. One last note is that our 
calculations do not include frequencies that have been observed only one time (as in PG 1219 and 
Feige 48) or stars for which we do not have the data. Both EC 14026 and EC 20338 are reported 
to have sufficient amplitude variability that frequencies completely disappear between observing 
seasons (Kilkenny et al. 2006b). 
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Figure 15. Mode density of pulsators organized by gravity. The dotted line indicates the limit of 3 frequencies per 1000 fiHz (all m = frequencies) 
and the dashed line indicates the limit of 9 frequencies per 1000 /jHz (all possible 21 + 1 m values) assuming £ ^ 2. 



Comparison with theoretical instability contours: We can make a direct comparison be 



tween the instability zone of the second-generation pulsation models of lCharpinet et al. 



(2001D and 



observed pulsation properties. Figures [FT] and [T8] are used for these discussions. In all the plots, 
filled circles are the sdBV stars of this study (stars indicate EC 20338, PG 0048, and PG 01540), 
open circles are sdB s tars (and inferred to be non-pulsators) from the Hamburg-Schmidt survey 



(Edelmann et al 



20031), Moehler et al (1990), and Saffer et al. (1994), and filled (blue) triangles 



2 We have inferred log g from the shortest pulsation frequency and T c ff from B — V and J — K s colours for EC 20338, PG 0048, and PG 0154. 
These are very crude estimates and should be considered as such. 
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Table 7. Mean amplitude, standard deviation, and their ratio for readily resolvable frequencies. Columns 6 through 8 provide the maximum and 
minimum measured amplitudes and the longest time covered by the observations. Except for HS 2201, all amplitudes and standard deviations have 
been calculated specifically for this paper from data we have. HS 2201 amplitudes are directly from Silvotti et al. 2002a and no error estimates 
were provided. Note that all frequencies are in fiHz and all amplitudes are in mma. 
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are PG 1716-type pulsators (Green et al. 2003). The contours are reproduced from CFB01 with the 
outside contour representing one unstable £ = frequency with each interior contour representing 
an additional unstable £ = frequency up to N = 7. 

In the left panel of Fig. [FTl we determine how discriminating the instability zone is by exam- 
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Figure 16. a(A)/(A) of individual frequencies. Triangles (squares) indicate phase-stable (-not-stable) frequencies, circles indicate frequencies 
with no phase information or indeterminate phase stability and stars indicate the PG 0048 frequencies which appear stochastic in nature. The 
dashed line is the JCD01 stochastic value. 



Table 8. Table 7 continued 
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PG1618B 7191 2.13 0.44 0.21 ±0.12 3.6 1.14 45 days 

7755 1.75 0.50 0.28 ±0.14 3.50 0.76 45 days 
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Figure 17. Comparison of observed sdB stars to pulsational instability contours (adapted from CFB01). Open circles are non-pulsating sdB stars, 
filled circles are EC 14026-type pulsators, and (blue) triangles are PG 1716-type pulsators. Left panel includes a dashed line indicating the inferred 
separation between EC 14026 and PG 1716-type pulsators and in the right panel the dot size corresponds to the number of total pulsation frequencies 
detected. 

ining the ratio of pulsators to non-pulsators within each contour. However, we need to have an 
exception for the red (cool) edge of the instability zone. There is an indication that sdB stars may 
switch from EC 14026- to PG 1716-type pulsators in this region. As such, this region (separated 
with a dotted line) should be excluded; particularly since it is stated that "all cool sdB stars of 
low gravity may be PG 1716 pulsators" (Fontaine et al. 2006). Working from the inside (N = 7) 
contour outward, the fraction of pulsating to non-pulsating sdB stars is: N ^ 7, 50%; N ^ 6, 
25%; N ^ 5, 25%; iV ^ 4, 22%; N ^ 3, 21%; N ^ 2, 20%; N ^ 1, 20% and the fraction 
of all sdB stars within the contours is iV ^ 7, 12%; N ^ Q, 52%; iV ^ 5, 80%, iV ^ 4, 90%; 
N ^ 3, 92%; N ^ 2, 95%; N ^ 1, 95%; N < 1, 100%. There does appear to be a relationship 
between the interior instability (N = 7) contour and fraction of pulsators. Yet outside of the first 
contour, the ratio only changes by 5% and all of the EC 14026-type pulsators are within the first 
three contours. Yet so are 80% of all sdB stars. So while the contours match where the pulsators 
are, they also match where most sdB stars with T cS ^30 000 K are. 

The remaining panels of Figs. [17] and [T8] show the same points, but the size of the dots cor- 
relate with the following observed properties; the number of pulsation frequencies (right panel 
of Fig. [171) . the maximum amplitude (A max ), and total pulsation power. A max and total pulsa- 
tion power are related as power is dominated by a few high-amplitude frequencies. Stars such 
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Figure 18. Same as in Fig.[T7]except that dot sizes for pulsators represent A max (left panel) and total pulsation power (right panel). 

as PG 0048 with many, but low amplitude frequencies cannot match the power of a single A ^ 
20 mma frequency. We would expect a correlation of the instability contours with the number 
of frequencies detected, as each interior contour increases the number of theoretically unstable 
frequencies, but such is not the case. From §5.2, we know that it is not a detection issue as rich 
pulsators occur with both low and high amplitudes. In fact there seems to be no correlation whatso- 
ever with the largest points (largest number of frequencies, highest amplitude and most pulsation 
power) occupying multiple regions of the log g — T eff diagrams and similar results for the smallest 
points. As such the group properties do not add observational support to the driving theory. 



6 CONCLUSIONS 

From extensive follow-up data acquired at MDM observatory, we are confident that we have re- 
solved the pulsation spectra of two additional pulsating sdB stars. For HS 0039, we detect 10 
additional frequencies bringing the total to 14 and for HS 0444, we confirm the two frequencies of 
the discovery data and detect an additional low amplitude frequency. We have also noted that while 
the amplitudes and phases of HS 0444 appear very steady over the duration of our observations, 
those in HS 0039 did not, but rather have varied considerably. This is illustrative of the variety ob- 
served in sdBV stars where some stars can have very simple and/or stable pulsation spectra while 
others can be quite rich, with tens of frequencies that may change amplitudes on short time scales. 
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Since the discovery of sdBV stars in 1996, more than 23 of the 34 known EC 14026-type pul- 
sators have received follow-up observations. We have examined these 23 stars for which extended 
timebase (and often multisite) observations have been acquired. We have searched for trends in 
frequency groupings, the number of high (H) amplitude to total (T) frequencies, and frequency 
density as a function of gravity and note that the only trend seems to be a weak relationship be- 
tween the H/T ratio and gravity: The lowest gravity pulsators have the smallest H/T ratio (a few 
very high amplitude frequencies) while the highest gravity pulsators have H/T^ 1 (relatively even 
amplitudes distributed amongst the observed frequencies). 

We examined amplitude stability, which has been used to infer stochastically excited oscil- 
lations in other variable classes and previously applied to sdB pulsators. And while many pul- 
sation frequencies fit the JCD01 value, the distribution between those and amplitude- stable and 
presumably driven is relatively smooth with no clear separations. However, phase-stable frequen- 
cies, some with quite large variability, preclude them from being stochastically excited and so the 
simplest conclusion is that the JCD01 criteria is not appropriate for sdBV stars. 

We compared pulsators to theoretical instability contours to search for relationships that cor- 
relate with current driving theory. While there is a weak concentration of pulsators precisely 
where expected, these same contours include more than 80% of all sdB stars in our sample with 
T eff ^ 30 000 K (an inferred cut-off for PG 1716-type pulsators) and there is no correlation be- 
tween the contours and the number of pulsation frequencies, nor their maximum amplitude or total 
pulsation power. It would be nice if we could note that any correlation is canceled by stars with 
lower gravities requiring less energy to drive higher amplitudes which results in more frequencies 
being detected, but this is also not the case as the amplitudes of neighboring stars anywhere in the 
log g — T eS plane can have amplitudes that differ by more than an order of magnitude. 

In the end, what we have uncovered in our examination of sdBV stars is the large variety they 
encompass. Subdwarf B pulsators seem to exist at all temperatures and gravities (though we only 
examined the EC 14026-type pulsators at T eff ^ 30 000 K) that sdB stars do and at roughly the 
same concentrations. Regardless of temperature or gravity (and thus evolutionary age and/or total 
mass and/or envelope mass and/or metallicity) pulsators can have relatively high (> 20 mma) or 
low (^2 mma), stable or unstable amplitudes and the frequency density can be quite high or very 
low. We hope that our investigation of these 23 well-studied sdBV stars will be useful and bring 
new insight to their theoretical perspectives. Additionally, we also hope that this study will be 
useful for observers pursuing multicolour photometry and time-series spectroscopy. We anticipate 
that such works will be required for an overall understanding of sdB pulsations, yet these more 
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detailed observations require that the pulsation characteristics of the star be known beforehand. So 
even though 23 of 34 known sdBV stars have received follow-up observations, there is still a long 
way to go both observationally and theoretically. 
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